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Abstract 



We show that the spectrum and eigenstates of open-shell multicharged 
atomic ions near the ionization threshold are chaotic, as a result of extremely 
high level densities of multiply excited electron states (10 3 eV _1 in Au 24+ ) and 
a strong configuration mixing. This complexity enables one to use statistical 
methods to analyze the system. The orbital occupation numbers obey the 
Fermi-Dirac distribution, and temperature can be introduced. We show that 
radiative capture of electrons through the multielectron resonant states is 
strongly enhanced compared to the direct radiative recombination. 

PACS numbers: 31.50,+w, 34.80.Lx, 32.70.Cs, 05.30.Fk 
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In this paper we investigate the spectrum and eigenstates of a multicharged positive ion 
at energies close to its ionization threshold /. Using Au 24+ (/ = 750 eV) as an example, 
we show that this spectrum is dominated by multiple electron excitations into a few low- 
lying unoccupied orbitals. As a result, it is extremely dense, with level spacings ~ 1 meV 
between the states of a given total angular momentum and parity J n . The electron Coulomb 
interaction induces a strong mixing of the multiply-excited configurations, which leads to 
a statistical equilibrium in the system. The latter is similar to a thermal equilibrium, and 
variables such as temperature can be introduced to describe it. This enables one to use 
statistical methods in the situation where a full dynamical quantum calculation is simply 
impossible because of the enormous size of the Hilbert space (•> 10 5 for Au 24+ ). 

We apply this approach to the problem of radiative capture of low-energy electrons 
by multicharged positive ions, and show that the contribution of resonant multielectronic 
recombination, which proceeds via electron capture into the multiply-excited compound 
states, is much greater than that of the direct radiative recombination. Our calculation gives 
a quantitative explanation of huge enhancements of the recombination rates, and removes 
the "enormous discrepancies between theoretical and experimental rate coefficients" The 
situation here turns out to be similar to the radiative neutron capture by complex nuclei 
[(n, 7) reaction] where the resonance mechanism involving the compound nucleus states is 
also much stronger than the direct capture 0. 

So far the enhanced recombination at low electron energies < 1 eV has been observed 
for a number of ions Q. Its magnitude ranges from a factor of about ten for Ar 13+ , Au 50+ , 
Pb 53+ , and U 28+ ||, to over a hundred for Au 25+ [p]]. This enhancement is sensitive to the 
electronic structure of the target, e.g., the recombination rates of Au 49+ and Au 51+ are much 
smaller than that of Au 50+ 0]. For few-electron ions, e.g., C 4+ , Ne 7+ and Ar 15+ [§-0], the 
observed rates are described well by the sum of the direct and dielectronic recombination 
rates. In more complicated cases, like U 28+ or Au 25+ , the questions of what are the particular 
resonances just above the threshold and how they contribute to the recombination "remain 
a mystery" ||. 

Let us consider the problem of electron recombination on Au 25+ . Due to electron cor- 
relations the electron can be captured into an excited state of the compound Au 24+ ion. 
Au 24+ has 55 electrons. Its ground state belongs to the Is 2 . . . 4/ 9 configuration. Figure [l] 
shows the energies of its relativistic orbitals nlj from a self-consistent Dirac-Fock calcula- 
tion. The energy of the highest orbital occupied in the ground state is £4/ 7/2 = —29.7 a.u. 
Our relativistic configuration-interaction (CI) calculation of the ground states of Au 24+ 4/ 9 
and Au 24+ 4/ 8 shows that they are characterized by J = y and 6, and their total energies 
are —18792.36 and —18764.80 a.u., respectively. Thus, the ionization threshold of Au 24+ is 
/ = 27.56 a.u.= 750 eV. 

The excited states of the ion are generated by transferring electrons from the ground state 
into the unoccupied orbitals above the Fermi level, or into the partially filled 4/ orbital. We 
are interested in the excitation spectrum of Au 24+ near its ionization threshold. This energy 
(27.5 a.u.) is sufficient to push up a few of the 4/ electrons, and even excite one or two 
electrons from the 4d orbital ||. Thus, we consider Au 24+ as a system of n = 19 electrons 
above the inactive Kr-like Is 2 . . . 4p 6 core. 

The number of multielectron states obtained by distributing 19 electrons over 31 rela- 
tivistic orbitals, 4<i 3/ / 2 through to 7gg/2, is enormous, even if we are only interested in the 
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excitation energies below 27.5 a.u. It is impossible to perform any CI calculation for them. 
However, there is another simpler way to analyze the spectrum. The scale of the configu- 
ration interaction strength is determined by the two-body Coulomb matrix elements. Their 
typical size in Au 24+ is about 1 a.u., i.e., much smaller than I. Configuration mixing aside, 
the CI does not shift the mean energies of the configurations. Therefore, we can construct 
the excitation spectrum of Au 24+ by calculating the mean energies Ei of the configurations, 

Ei = E core + 2^ e a n a + ^ — — — U ab: (1) 

a<b 1 + °ab 

and the numbers of many-electron states in each of them, iVj = Y[ a gJ/[nJ(g a ~ n aV-}, where 
n a are the orbital occupation numbers in the configuration (S a n a = n), e a is the single- 
particle energy of the orbital a in the field of the core, g a = 2j a + 1, and U a b is the average 
Coulomb interaction between the electrons in orbitals a and b (direct — exchange). 

We find that there are 9000 configurations within 35 a.u. of the Au 24+ ground state. 
They comprise a total of 2.1 x 10 8 many-electron states, which corresponds to about 5 x 10 5 
states in each J 77 manifold. Figure |2| shows the total density of states 

p(E)=J2N l 6(E-E l ) (2) 



averaged over 1 a.u. energy intervals, as a function of v E, where E is the excitation energy 
above the ground state. The inset presents a break-up of the total density near the ionization 
threshold into the densities of states with given J: p(E) = 2j(2J + l)pj(E). The most 
abundant values are J = § _ y- For a given parity the density of such states at E ~ I is 
Pjk « 3.5 x 10 4 a.u., which corresponds to a mean level spacing D = 1/pj* ~ 1 meV. Figure 
H demonstrates the characteristic p oc exp(ay/E) behaviour of the level density predicted by 
the Fermi-gas model plj , where a is related to the single-particle level density at the Fermi 
level g{eF) = 3a 2 /2ir 2 . A Fermi-gas model ansatz 

p(E) = AE~ V exp(av / ^) , (3) 

with A = 31.6, v = 1.56, and a = 3.35 gives an accurate fit of the level density at E > 1 
a.u. The corresponding g(ep) = 1.7 a.u. is in agreement with the orbital spectrum in Fig. 
[j]. For most abundant J n states pjn(E) is given by Eq. (^|) with Ajn w 0.15. 

At first sight the huge level density makes the Au 24+ problem very complicated. In 
reality this complexity enables one to use statistical methods to describe the system. The 
interaction between multiply-excited configuration states mixes them completely, and they 
loose their individual features. In this regime the spectral statistics are close to those 
of a random matrix ensemble, the eigenstates cannot be characterized by any quantum 
numbers except the exact ones (energy and J n ), and the orbital occupation numbers deviate 
prominently from integers. This regime can be described as many-body quantum chaos. We 
have extensively studied it in direct numerical calculations for the rare-earth atom of Ce - 



a system with four valence electrons 11-13 



The strength of the configuration mixing is characterized by the spreading width r spr . 
For a configuration basis state $fc with energy E^ it defines the energy range \E — E^\ ^ r spr 
of eigenstates in which this basis state noticeably participates. By the same token it shows 
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that a particular eigenstate \1> = J2kCk$k contains a large number N ~ T spT /D of principal 
components - basis states characterized by Ck ~ 1/ y~N. Outside the spreading width Ck 



decrease. Apart from this, Ck behave like random variables |TTJ. The effect of spreading is 
approximated well by the Breit-Wigner shape [T(J 



The normalization C\ = 1 yields iV = ttT/2D. In systems with small level spacings D the 
number of principal components N can be very large. It reaches several hundreds in Ce, and 
10 6 in complex nuclei. In Fig. § we illustrate this behaviour by the results of a CI calculation 
for two Au 24+ configurations near the ionization threshold: 4/|/ 2 4/^ 2 5pi/25p3/25/7/2 and 

4/5/24/|/25pi/25<i3/25g7/2- These two configurations produce 1254 J w = y + states. Their 
mixing is complete, since the weight of each configuration in every eigenstate is about 50%. 
A Breit-Wigner fit of the mean-squared components yields N = 975 and r spr = 0.50 a.u. 
The spreading width is related to the mean-squared off-diagonal Hamiltonian matrix element 
and the mean level spacing as r spr ~ 2ttH^/ D [[l 7 ]]. When more configurations are included, 

both D and Hfj decrease, and r spr does not change much. If one were able to do a full-scale 
CI calculation near the ionization threshold of Au 24+ eigenstates with N = (7r/2)T 8pT pj^ ~ 
3 x 10 4 would be obtained. 

The spreading of the basis states due to configuration interaction introduces a natural 
statistical averaging in the system. Based on this averaging, a statistical theory of finite 
Fermi systems of interacting particles can be developed [14]. It enables one to calculate 



various properties of the system without diagonalizing the Hamiltonian, as sums over the 
basis states, e.g., the mean occupations numbers, n a (E) = ^2 k Ck(E)n^' > where is the 
occupation number of the orbital a in the basis state k. Using a simple Gaussian model 
spreading we calculate mean orbital occupation numbers at different excitation energies. 
Figure |4] shows that the result (circles) is described well by the Fermi-Dirac formula n a = 
{1 + exp[(e a — /x)/T]} _1 . The temperature T and the chemical potential p are chosen to 
give the best fit of the numerical occupation numbers. On the other hand, the temperature 
can be determined from the level density p{E), Eq. (0), through the canonical average 
E{T) = Z^ 1 J e~ E l 7 'Ep(E)dE, where Z = J e~ E l T p(E)dE, or from a statistical formula 
T _1 = d\n[p(E)]/dE, using the smooth fit (|3]). The latter yields T ~ 2^/E/a, characteristic 
of Fermi systems. For energies above 3 a.u. all three definitions give close values, Fig. |5]. 

Let us now estimate the direct and resonant contributions to the recombination rate 
of Au 25+ . The direct radiative recombination cross section is estimated by introducing an 



effective ionic charge Z{ into the formula of Bethe and Salpeter |15 

where e is the electron energy and n is the principal quantum number of the lowest un- 
occupied ionic orbital (we use atomic units). Using Z 4 = 25, n = 5, and e = 0.1 eV we 
obtain erW »s 7 x 10~ 17 cm 2 , which gives a rate of A = av = 1.3 x 10 -9 cm 3 s _1 , two orders 
of magnitude smaller than the experimental A = 1.8 x 10~ 7 cm 3 s~ 1 at this energy p]]. 
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The energy-averaged cross section of the resonant radiative capture of a low-energy s 
electron is IlIBI 



7T 2 r r 7T 2 r 

e D(r 7 + r e ) £ d ^ ie>>i ^' 

where T 7 and r e are the mean radiative and autoionization (or elastic) widths of the reso- 
nances, D is the mean resonance spacing, and we drop the statistical weights of the initial 
and intermediate ionic states. The relation T e ^> T 7 is usually valid for a few lower partial 
waves, where the electron interaction is stronger than the electromagnetic one. 

The radiative width of the resonant state at energy E w / is found by summing the 
partial widths for all lower- lying states E' = E — u, 

r ^27TT7o -£^pMi -<*)*<>, (7) 

where the factor 3 accounts for J 1 = J, J ± 1, and d w is the reduced dipole matrix element 
between the many-electron states. Because of the chaotic structure of these states is 
suppressed compared to the typical single-particle matrix element do- d w ~ d /\/N P, [lT| , [l2 



This estimate for systems with dense chaotic spectra in fact follows from the dipole sum rule: 
the number of lines in the spectrum is large, oc D^ 1 oc N, consequently, the line strengths 
are small, \d u \ 2 ~ {d^N^ 1 . 

Using this estimate and calculating the integral in Eq. (0) by the saddle-point method 
we obtain 

where lu = Q\^I/a corresponds to the maximum of the the decay photon spectrum in Eq. 
(0). This cross section has the same energy dependence as a^ d \ To evaluate its magnitude 
we use do ~ , 2J+1 rs 10, and substitute r spr = 0.5, too = 9.4, and pj^{I—ujq) = 2.5 xlO 3 
a.u. At e = 0.1 eV this gives a^' = 7 x 10 -16 cm 2 , therefore, /a^ A > = 10, and we obtain 
a factor of ten enhancement over the direct recombination. It comes from the large effective 
number of final states in the radiative width in Eq. (^) (numerically T 7 ~ 2 x 10~ 7 a.u.). 
If we include contributions of the higher electron partial waves the enhancement will match 
the experimentally observed values. 

In summary, the resonant radiative capture mechanism fully explains the strongly en- 
hanced recombination rates observed for eV electrons on multicharged ions. Its origin is in 
the high level densities of chaotic multiply-excited electron states in multicharged ions. The 
size of the enhancement is sensitive to the electron structure of the ion, which determines 
the level density. We have shown that a statistical approach can be applied to the analysis 
of this complex system. One can also use a statistical theory to calculate mean-squared ma- 
trix elements between multiply excited chaotic states in terms of single-particle amplitudes, 
occupation numbers, r spr and D JTT| , p!2] , pT| , and obtain accurate quantitative information 



about the processes involving chaotic states and resonances. At higher electron energies the 
resonant capture proceeds via so-called doorway states [ID] - simple dielectronic autoionizing 



states, which are then "fragmented" into the dense spectrum of multiply-excited resonances 
(see |§ and JE2| and Refs. therein). 



5 



REFERENCES 



[1] A. Hoffknecht et ai, J. Phys. B 31, 2415 (1998). 

[2] V. V. Flambaum and 0. P. Sushkov, Nuclear Physics A412, 13 (1984); A435, 352 
(1985). 

[3] We do not consider a specific enhancement below 1 meV observed for all ions including 
the fully stripped Ar 18+ , O. Uwira et ai, Hyperfine Interact. 108, 167 (1997). 

[4] H. Gao et ai, Phys. Rev. Lett. 75, 4381 (1995); O. Uwira et ai, Hyperfine Interact. 
108, 149 (1997); A. Miiller and A. Wolf, ibid. 109, 233 (1997). 

[5] S. Schennach et ai, Z. Phys. D 30, 291 (1994). 

[6] W. Zong et al, Phys. Rev. A 56, 386 (1997). 

[7] R. Schuch et al, Phys. Scripta T73, 114 (1997). 

[8] D. M. Mitnik et al, Phys. Rev. A 57, 4365 (1998). 

[9] Infinite Rydberg series belong to the single-particle aspect of the problem, and we do 
not consider them here. 
[10] A. Bohr and B. Mottelson, Nuclear structure, Vol. 1 (Benjamin, New York, 1969). 
[11] V. V. Flambaum et al, Phys. Rev. A 50, 267 (1994). 

[12] V. V. Flambaum, A. A. Gribakina, and G. F. Gribakin, Phys. Rev. A 54, 2066 (1996); 

58, 230 (1998). 
[13] V. V. Flambaum et ai, Phys. Rev. E 57, 4933 (1998). 

[14] V. V. Flambaum and F. M. Izrailev, Phys. Rev. E 55 R13 (1997); 56 5144 (1997). 
[15] H. Bethe and E. Salpeter, Quantum Mechanics of Atoms with One and Two Electrons 

(Springer, Berlin, 1957). 
[16] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, (Pergamon Press, Oxford, 1977), 

Ch. 18. 

[17] V. V. Flambaum and O. K. Vorov, Phys. Rev. Lett. 70, 4051 (1993). 



6 



FIGURES 



10 



20 



too 

CD 

W-40 



50 



60 




4s 







4f 



-1c 



4p 



J I L 

1 2 3 

Orbital momentum 



4 



FIG. 1. Electron orbitals of Au 24+ from the Dirac-Fock calculation of its ground state 
ls 2 ...4/ 9 . 
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FIG. 3. Components of a 590th J n = eigenstate from a two-configuration calculation (top), 
and a fit of C%(E) by the Breit-Wigner formula (bottom). 
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FIG. 4. Orbital occupation numbers in Au 24+ calculated at E = 1, 9.5, and 27.5 a.u. (circles), 
compared with the Fermi-Dirac distribution (solid line). 
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FIG. 5. Temperature vs energy for Au 24+ . Solid line - canonical definition; dashed line 
statistical definition; solid circles - Fermi-Dirac fits of the occupation numbers. 
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